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*5 : ABSTRACT 

6 : 

We have identified 63 flares on M dwarfs from the individual component 
spectra in the Sloan Digital Sky Survey using a novel measurement of emission 
line strength called the Flare Line Index. Each of the ~38,000 M dwarfs in the 
SDSS low mass star spectroscopic sample of West et al. was observed several 
times (usually 3-5) in exposures that were typically 9-25 minutes in duration. 
Our criteria allowed us to identify flares that exhibit very strong Ha and H/3 
i ; emission line strength and / or significant variability in those lines throughout the 

course of the exposures. The flares we identified have characteristics consistent 
with flares observed by classical spectroscopic monitoring. The flare duty cycle 
for the objects in our sample is found to increase from 0.02% for early M dwarfs to 
. ^ ; 3% for late M dwarfs. We find that the flare duty cycle is larger in the population 

near the Galactic plane and that the flare stars are more spatially restricted than 
the magnetically active but non-flaring stars. This suggests that flare frequency 
may be related to stellar age (younger stars are more likely to flare) and that the 
flare stars are younger than the mean active population. 

Subject headings: stars: flare, stars: late-type 



1. Introduction 



Flares are exp losive events caused by magnetic reconnection in stellar atmospheres 
( IHaisch et al.lll99ll and references therein). In the standard model, electrons are accelerated 
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along magnetic field lines and impact the lower stellar atmosphere. This catastrophic r elease 



of magnetic energy causes emission from t he radio to the X- ray (IHawley et al.l Il995[ 12003 



Osten et al.ll2005t iFuhrmeister et al.l 120071 ; iBerger et al.ll2008l ). On M dwarfs, the most dis- 
tinct observational characteristic is the tremendous increase in the blue and near-ultraviolet 



continuum emission, up t o several magni t udes in a few minutes or less (IHawley fc Pettersen 
199ll ; lEason et al.l Il992l ; iKowalski et al.l 120101 ) . The initial burst, or impulsive phase, is 
also characterized in the optical spectru m by increases in chromospheric line emission, par 



Fuhrmeister fc Schmitt 



ticularly in the hydrogen Ba lmer lines (IHawley fc Pettersenl Il99ll ; iMartm fc Ardilal 12001 



20041 ). The impulsive phase is followe d by a gradual decay phase 



which can last from tens of minutes to hours for the largest flares (lMoffettlll974l ; iZhilyaev et al. 
200/1 



The rate at which flares occur is of significant current interest. Recently, numer- 
ous pla nets have been d i scovered around M dwarfs, including some o f only a few Earth 



masses (lUdrv et all 120071: iCharbonneau et al.ll2009l ; iForveille et al.ll2009l ; iMavor et al.l 12009 
Correia et al. 2Qloh . The atmospheres of these planets may be greatly affected by the amount 



of high-energy radiation incident upon them. Cool, red M dwarfs produce little high-energy 
radiation du ring quiescence, so the flaring rate is an important facto r in determining planet 



Segura et al 



201 



habitabilitv (iHeath et al.lll999l ; ISegura et al.ll2005l ; iTarter et al.ll2007l ; IWalkowicz et al 



2008 



Additionally, char acterizing th e M dwarf flaring rate is impor tant for new ins truments 
such as Pan-STARRS (Kaiser 2004), the Palomar Transit Factory jRau et al. 20091) and th e 
Large Synoptic Survey Telescope ( IlSST Science Collaborations: Paul A. Abell et al. 2009 ). 
which will carry out large, all-sky surveys in the time-domain. Effectively selecting rare, 
exotic transients from the large sample of M dwarf flares that will be observed in these 
surveys requires reliable knowledge of expected flare characteristics and rates. 

Traditionally, the frequency of flares on M dw arfs has been det ermi ned through time- 

resolv ed photometric monitoring of individual stars. lLacy et al.l (119761 ) and lGershberg fe Shakhovskaia 
(119831 ) have amassed several hundred hours of optical observations on dozens of the most 
magnetically active and well-known flare stars and have derived power law relationships be- 
tween the total energy in the U band and the frequency of a flare, with less energetic flares 
occurring more frequently. Studies on individual stars have co nfirmed the power law form of 
the f r equency distribution, albeit with a range of exponents 



1984; Robinson et al. 



995 



Leto et al. 



1997; 



Robinson et al. 



lengths, lAudard et al.l (120001 ) and ISanz-Forcada & Micelal (120021 ) found that t he EUVE flare 



Walker 1981; Pettersen et al. 



9991 ). At ultraviolet wave- 



frequency distributions for active stars also had a power law form. In addition, lAudard et al. 



(120001 ) demonstrated that large flares occur preferentially on the X-ray brightest stars. 



-3 - 



Spectroscopic studies of flares have also investigated the flare rate, but with smaller 
samples and incomplete information on flare energies it is not possible to compute the flare 
frequency distribution as described above for photometric data. Instead the results are 
usually presented as "flare duty cycles" , meaning that during the entire period of observation, 
flares occur during a certain percentage of the time. Alternatively a "flare rate" (number 
of flares per hour) is sometimes given, but without reference to the energy distribution of 
the flares. Both the flare duty cycle and flare rate are strong functions of the observational 
sample, since such factors as the individual exposure time (compared to the typical flare 
timescale), the flare visibility in the line radiation compared to the continuum, and the flare 
energy will all be important in determining whether a particular observation is counted as 
flaring. 



Reid et al.l (119991 ) estimated a duty cycl e of 7% from Ha observations of the M9.5 dwarf 
BRI 0021-0214, while ICrespo-Chacon et al.l (120061 ) observed 14 small flares on AD Leonis, 
none of which produced measurable continuum increase, and determined a flare rate of 



> 0.7 flares/hr. ISchmidt et al.l (120071 ) investigated a sample of 81 M dwarfs, each observed a 
few ti mes and found a flare duty cycle of ~5%, based on Ha variability. Recently, iLee et al. 



(120101 ) reported a duty cycle of < 5% for Ha flares (defined as a brightness increase of a 
factor of ten) from time-resolved spectroscopic monitoring (~1 hour each) of 43 M dwarfs. 

The small numbers of stars in most previous studies (both photometric and spectro- 
scopic) did not allow e xamination of whether the flare rate is a function of stellar or Galactic 
parameters. However, iKowalski et al.l ( 120091 ) used the repeat observations from the SDSS 
Stripe 82 photometry (~60 epochs per object spread over 10 years) to identify flares and 
determine flare rates. They identified 271 flares on more than 50,000 M dwarfs, a sample 
which is orders of magnitude larger than previous work and the first to contain significant 
numbers o f both magnetically ac tive and inactive stars. The vast majority of flares (but not 
all) in the IKowalski et al.l ( 120091 ) study occurred on active stars, and they found that the 
flare star spatial distribut ions reflect the active star distributions as a function of spectral 
subtype rtWest et al.ll2008h . 



It should be clear that the flare rates, flare duty cycles, and flare frequency distribu- 
tions described above are distinctly different measurements, and are difficult to compare. In 
addition, the criteria used to define a flare differ depending upon the method of observation 
(photometry vs spectroscopy), the exposure time during which a flare may occur (resolved 
vs unresolved), the wavelength (optical vs near-ultraviolet), and the emission mechanism 
(continuum vs lines). The sample selection (spectral type distribution, Galactic height dis- 
tribution, color selection, etc.) also has a strong influence on the number of flares found. 
Additionally, flares themselves are not homogeneous. Flares with the same energy exhibit 



-4- 



apparent magnitude variations that depend on the brightness of the star during quiescence 
(the so-called "contrast" effect). On stars with the same quiescent brightness, flares of the 
same energy exhibit a wide range of light curve shapes, with some flares rising and decaying 
quickly, giving large changes in apparent magnitude, while others are slower, with smaller 
peak magnitude. Understanding how often an M dwarf is seen in a flaring state requires 
knowledge of many stellar, flare, and sample properties. 

Spectroscopic flare data typically have longer exposure times, may contain unresolved 
flares, are usually at optical wavelengths, and are more likely to observe flares, especially 
small ones, through enhanced line emission. Flare rates and frequency distributions are 
therefore not good descriptions of the data, and we will use the flare duty cycle (the fraction 
of epochs classified as flares) in this work. 

Here we present the first large, statistical study of spectroscopic flare duty cycles. The 
low mass st ars in this study are taken from the Sloan Digital Sky Survey Dat a Release 5 
(SPS S DR5: lAdelman-McCarthv et al.ll2007h low-mass star spectroscopic sample (jWest et al.l 



20081 ). which consists of over 38,000 objects between spectral types M0 and L0. All SDSS 
spectra are the result of coadding several (usually 3-5) individual exposures, typically (but 
not always) observed on the same night. The individual component exposure times were 
generally between 9 a nd 25 minutes. The individual sp ectra were made public as part of the 



SDSS Data Release 6 (lAdelman-McCarthy et al.ll2008l ). The repeated, shorter exposures are 



used to aid in cosmic ray rejection for the final combined SDSS spectrum. However, these 
repeat spectra can be used to examine the time variability of spectral features. We take 
advantage of the time sampling contained in the individual exposures to identify flares based 
on the strength and variability of spectral lines. 

We have created a unique tool we use to identify flares: the Flare Line Index (FLI), 
which allows us to quantitatively measure spectral line strength in low signal-to-noise (SNR) 
spectra, where eq uivalent width mea surements are not possible. We then used False Discov- 



ery Rate analysis (IMiller et al.ll200lf ) to identify flares based on exceptionally strong emission 
line strength. Additionally, we used FLI variability criteria to identify flares based on changes 
in emission line strength over the course of the exposures. 

In §|2] we further characterize our sample and describe our quality cuts. The tools for 
our analysis are described in detail in £j3j The flare selection criteria are defined in §HJ and 
the results are presented in £j5j Finally, we conclude in §|6j 
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2. The Sample 



We began with the SDSS DR5 low-mass star sample of IWest et al.l (120081 ). for which 
low mass stellar candidates were chosen from the SDSS DR5 spectroscopic catalog using 
the photometric color selection r — i > 0.53 and % — z > 0.3, along with a standard set 
of SDSS processing flags (SATURATED, BRIGHT, NODEBLEND, INTERP_CENTER, 
BAD_COUNTS_ERROR, PEAKCENTER, NOTCHECKED, and NOPROFILE were all re- 
quired to be zero). Stars with r-band extinction > 0.5, colors consistent with having a white 
dwarf companion (ISmolcic et al.1 2004'). or measured radial velocity > 500 km s _1 were culled 



from the sample (see IWest et al 



( 120081 ) for more details). 



The 38,717 spectra of individual objects in the DR5 sample that meet the above criteria 
consist of 142,032 individual exposures that were made available in DR6. The distribution 
of number of exposures per object and the exposure times are shown in Figure [TJ Although 
most of the spectra for an object were observed consecutively on a single night, 27% of 
objects have additional exposures separated by two hours or more. 

The SDSS spectra have wavelength coverage from 3800 to 9200 A, and resolution ranging 
from 1800 to 2200. We therefore were able to measure five hydrogen Balmer lines (H8, H<5,H_7, 
H/3, Ha), as well as Ca II K. Ca II H and He are blended at this resolution and are not 
considered in our analy sis. The wavelength re gions used for the continuum and line regions 
are given in Table 1 (cf. iBochanski et al.ll2007l ). An example of the data, Figure [2] shows the 
four exposures and line profiles from the six measured lines of a typical decay phase flare. 



Table 1. Line and continuum regions 



Line 


Line Re 


gion (A) 


Cont 


inuum 1 (A) 


Cont 


inuum 2 (A) 


H8 


3884.15 


- 3898.15 


3850 


- 3880 


3910 


- 3930 


Ca II K 


3928.66 


- 3938.66 


3910 


- 3925 


3950 


- 3960 


E5 


4097.00 


- 4110.00 


4030 


- 4080 


4120 


- 4170 


H 7 


4331.69 


- 4350.00 


4270 


- 4320 


4360 


- 4410 


H/3 


4855.72 


- 4870.00 


4810 


- 4850 


4880 


- 4900 


Ha 


6557.61 


- 6571.61 


6500 


- 6550 


6575 


- 6625 
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2.1. Quality Cuts 

2.1.1. Cosmic Rays 

Cosmic Rays (CRs) were not removed from the individual component spectra in the 
SDSS DR6 spectroscopic pipeline, but individual pixels that contained likely CRs were 
flagged. We visually inspected many of the flagged pixels and found that the pipeline rou- 
tines frequently misidentified stellar emission lines as CRs. Since flares produce strong and 
time variable emission that might be misclassified as CRs, we developed a separate crite- 
rion to identify pixels as CRs in the line regions, so as not to reject genuine emission: a 
pixel within the line region was called a CR if the difference between that pixel and the 
continuum level was more than ten times the standard deviation of the other pixels in the 
region. Our criterion identified CRs in the Ha and H/3 line regions in 0.55% and 0.69% of 
spectra respectively, as compared to 2.7% and 3.5% of the spectra that the SDSS pipeline 
flags identified. There were no significant variations in the rate of CRs detected at different 
spectral subtypes using either definition. 

Although this criterion still allows some CRs to remain undetected and masquerade 
as emission lines, we require emission in both Ha and H/3 to identify an object as in a 
flaring state (see $3J). This requirement greatly reduces the contamination by CRs, since 
two independent CR strikes on H lines in the same spectrum happen very rarely. Using 
the SDSS CR flags, we found that a "CR" is detected in both the Ha and H/3 regions in 
0.3% of the spectra, compared to 0.007% with our CR criterion. All of the objects that 
met our flare criteria were inspected by eye, ensuring that no CRs led to an incorrect flare 
classification. Since the continuum regions are much wider than the line regions and are 
only used as an indicator of the continuum level and signal-to- noise ratio (SNR), we used 
the more conservative SDSS CR flags to mask any suspect pixels in all calculations involving 
the continuum regions. 

2.1.2. Sample Separation by Signal-to-Noise Ratio 

The large sample size allowed us to divide our objects into 3 quality bins based on the 
SNR. Because M dwarfs are red and often have low SNR in the blue portions of their spectra, 
our division into three bins identified those spectra that could be analyzed in both the Ha 
and H/3 regions. These divisions were also important for our False Discovery Rate analysis 
(see p . 

We made our cuts based on the continuum level in regions near the lines of interest 
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because accurate emission line measurements depend on our ability to define the local con- 
tinuum. The high signal-to-noise (HSN) sample includes objects where all component spectra 
have median continuum levels at least three times the standard deviation of the points in 
the continuum regions for both Ha and H/3. There are 43,788 spectra in the HSN sample, 
from 12,459 individual stars. Because earlier-type stars are both bluer and brighter, the 
HSN objects are almost entirely spectral types MO-4. The spectral type distribution of HSN 
spectra is shown in Figure [3j 

Our medium signal-to-noise (MSN) sample is defined in a similar way to the HSN sample, 
except that we only require the median Ha continuum level to be at least 3 times the standard 
deviation of the points in the region. The MSN spectra have very noisy continuum flux near 
H/3. However, many of these objects exhibited H/3 line emission that can be used in our 
flare analysis. Like the HSN sample, all of the spectra for each object were required to meet 
the MSN criteria in order to be included in the MSN sample. There are 68,203 spectra 
from 18,512 individual stars in the MSN sample. The spectral type distribution of the MSN 
spectra is also shown in Figure [3j 

For completeness, we note that for ~20% of the stars in the original sample, the indi- 
vidual component spectra have low SNR that did not meet the MSN cut. This LSN sample 
may well contain flares but cannot be used in the detailed analysis presented here. 



2.2. Spectral Type Determination 



The spectral type of each star was determined by IWest et al.l (120081) fro m the co-added 
DR5 spectra using the Hammer spectral typing facility (jCovey et al.l 120071 ). The Hammer 
uses measurements of molecular bands and line strengths to estim ate a spectral type and is 
accurate to within one subtype. We adopted the I West et al.l (120081 ) types for our analysis of 
the individual component spectra, which have lower SNR than the co-added spectra. 

The Hammer is tuned to identify and classify M dwarfs during quiescence by fitting 
templates to molecular band depths and overall spectral shape. Large flares produce optical 
continuum emission that veils the molecular features and makes the spectral shape signifi- 
cantly bluer. Therefore, we verified that the sample does not exclude continuum flares that 
cause the star to be misidentified by the Hammer. We spectral type d the dM4.5e flare sta r 
YZ CMi during quiescence and during a giant flare of AU ~ -6 mags (IKowalski et al.ll2010l ). 
The quiescent spectral type returned by the Hammer was M4, while during some stages 
of the flare, types of M2 and M3 were found, still reasonably close to the quiet value. To 
further establish that the spectral type determination was not biased against flares, we also 
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searched for flares in all spectroscopic objects whose photo metric colors (taken at a different 



time than the spectroscopy) were consistent with M dwarfs ( jWest et al.ll2008t iKowalski et al. 



20091 ). but were not identified as M dwarfs by the Hammer. There were no large, continuum- 



enhancement flares identified among these objects. These two tests confirm that we did not 
select against continuum flares by relying on automatic spectral typing. 



3. Flare Analysis Tools 



The unique SDSS data set provides an opportunity to observe the spectroscopic time 
evolution of a large number of flares. It does, however, present some challenges. First, the 
total length of consecutive observations is generally short (~45 minutes), and if a flare does 
occur, it may start before or at any time during the individual observations. Second, the 
time resolution is coarse: the duration of small flares is often less than the exposure time 
of one observation. Both of these factors make determining the time of the flare beginning, 
peak, or end ambiguous. Additionally, only the largest flares, which are rare, will be evident 
throughout an entire sequence of spectra. In the case of large flares that last for several 
hours, the observations may only cover a small fraction of the entire flare light curve. Finally, 
flares are more likely to be observed during the longer gradual phase, when the continuum 
enhancement is likely to be significantly diminished. On the other hand, the emission lines 
will remain enhanced throughout the flare. Given these challenges, we chose to concentrate 
on identifying flares based on strong and/or variable chromospheric line emission. 

Ha ha s been studied extensively in low mass stars, and is well known to vary on short 
timescales (IBopp fc Schmitall978l ; iGizis et al.ll2002l ; ICincunegui et al.ll2007l ; IWalkowicz fc Hawley 



2009; Lee et al. 



2010). even outside of flares. Thus, if Ha is solely used to identify flares 



there is no way of distinguishing variability from flaring (except for the largest Ha flares), 
given the small number of exposures and the relatively long exposure times in our sample. 
Although H/3 variability has not been as well-studied on M dwarfs, any flare large enough to 
be seen in a 15-minute exposure shou ld also show an increas e in H/3 emission, as evidenced 
by studies of individual flares, e.g. ([Martin fc Ardila 2001 ; Fuhrmeister fc Schmitt 2004 ; 
Crespo-Chacon et al. 2006 ; Kowalski et al] 2010 ). Therefore we required both the Ha and 
H/3 lines to meet our flare criteria for line strength and/or variability (see §H]for details of 
the flare criteria). This requirement has the additional benefit of ensuring that cosmic rays 
are not falsely classified as flares. 
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3.1. Flare Line Index 

A commonly used measurement of emission line strength is the equivalent width (EW). 
However, EW measurements are problematic when the continuum flux near the line is weak 
and has low SNR. Figure H] shows the Ha and H/3 line regions during a sequence of exposures 
of an M5 dwarf from our MSN sample. The H/3 continuum flux is very close to zero, so 
division by the continuum to obtain EWs gives values of -3A, 106A, and 19A for the three 
consecutive exposures. While the star in Figure H] is clearly not in a flaring state, the EW 
measurements for H/3 are large and variable, as they might be during a flare. This issue 
could be avoided by limiting the sample to objects with high SNR in the H/3 continuum (i.e., 
the HSN sample). However, including the MSN sample greatly increases our sample size. 
We therefore developed a new method for measuring line strength that is independent of the 
continuum strength. We define the Flare Line Index (FLI) as: 

FLI = l/a, (1) 

where I is the mean value of the flux in the line region (for a specific emission line) minus 
the continuum flux (associated with the emission line; see Table 1) and a is the standard 
deviation of the continuum. The FLI measures the line strength in terms of its significance 
compared to the noise. 

The FLI allows us to determine whether a line measurement is statistically significant 
even when the continuum is weak and/or uncertain. For the star in Figure HJ the FLI values 
for H/3 are -0.3, 1.1, and 0.2, respectively. The differences in these values can be attributed 
to the noise in the spectra. The FLI values are neither as large nor as variable as the EW 
measurements of the same data. When the continuum is well-measured, as in the Ha line, 
the FLI values and EW measurements are similar. The FLI does depend on the SNR of the 
spectra; stars with higher SNR have larger FLI values for the same emission line flux. 

FLI values for both Ha and Hf3 were computed for all spectra in the HSN and MSN 
samples and used in the subsequent analysis. 

3.2. False Discovery Rate Analysis 

Determining if an object is in a flaring state from a single spectrum is difficult because 
active stars have a range of nominal, non-flaring levels of Ha activity. Also, there may be 
little difference between the emission line strengths of an active star in quiescence and those 
of a flaring star during a small flare or in the gradual phase of a large flare. We searched for 
flares based on the strength of the line emission, and on the variability in repeated spectral 
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observations. We did not require that both criteria be met, since small flares may not have 
particularly strong lines, while large flares may not show much variability over the timescale 
of the SDSS observations. 



We used False Discovery Rate (FDR; iMiller et al.l 120011 ) analysis to search for spectra 



with particularly strong emission lines that were outliers from the typical stellar distribution 
of FLI values. Our sample exhibits a continuous distribution of measured FLI values for 
both Ha and H/3, ranging from lines in absorption through emission and into the flare 
regime. FDR analysis is a quantitative method for determining a threshold value of a test 
statistic that separates real sources from a background (null) distribution for a given false- 
discovery rate. In our case, we sought to determine the FLI value above which stars were 
considered to be in a flaring state. The analysis uses an adaptive technique that establishes 
the threshold value based on the desired false positive rate, defined by the parameter a. Our 
FDR method replaces an arbitrarily chosen confidence level (such as 3<r) with a known false 
positive rate and can select the outliers from our sample (as compared to the null) as flares 
while controlling the rate of false detections through the choice of a. For our analysis, we 
limited the false positive rate to a = 10%. 

The key to using FDR analysis is to assemble a meaningful null distribution. In the 
case of identifying flares, the null distribution should be a set of spectra that have similar 
properties to our candidate flare sample, but contain no flares. Active M dwarf spectra (that 
are not flaring) will provide a good null sample for the flare analysis. 



3.3. Constructing an Active Sample using FDR 



The fraction of stars that are ma gnetically active increases from a few percent at spectral 



v_p 

type MO to nearly 100% at type M9 flHawlev et al.lll996l : iGizis et al.ll200~of: Iwest et al.ll2004 



20081 ). Rather than simply adopting the active sample found by I West et al.l (120081 ). we tested 
the effectiveness of the FLI values and F DR analysis t ools o n our HSN and MSN samples 
by verifying that we could reproduce the I West et al.l (120081 ) results. Again, we require a 
null distribution to separate the active stars from the inac tive ones. For th is case, we 
adopted the DR5 composite spectra of inactive stars found by I West et al.l (120081 ) as the null 
distribution. We used the individual component spectra in the HSN and MSN samples as 
the test distribution. We carried out the analysis separately for each spectral type, although 
we found that the HSN sample only had sufficient numbers for the MO-3 types and the MSN 
sample for the MO-7 types (see Figure [3]). 



An example of the FDR analysis for M3 stars is shown in Figure The gray line 
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represents the distribution of FLI values for the Ha emission of the M3 dwarfs in the HSN 
sample (the test sample, which includes both active and inactive stars), while the black line 
represents the FLI distribution of the inactive M3 dwarfs from the DR5 composite spectra 
(the null sample, which contains only inactive stars). The distributions have been normalized 
to have the same peak value centered at FLI=0. The vertical line (at a FLI of 0.66) is the 
threshold value determined by the FDR analysis that divides active stars from inactive ones, 
with the false positive rate a set to 10%. The FLI emission threshold values we determined 
for the MO-3 spectral types in the HSN sample are shown in Figure [6] (labeled "Emission"). 
Note that for later types, the M3 value was adopted. Table 2 gives the Ha FLI emission 
threshold values for both the HSN and MSN samples. 

Using our FLI/FDR determination of activity, we found that 2.8% of the MOs are active 
and this fr action rises t o 84% of M7-9 stars. These compare well to the active fractions 



reported by lWest et al.l (120081 ). who found 2.6% and 90%, respectively. Although this agree- 
ment might be expected because we used the DR5 inactive spectra as our null distribution, 
we used a different test statistic (FLI instead of EW) and found our active sample using 
FDR analysis instead of defining an E W activity thresh old of 1A on statistically significant 



emission lines in high SNR spectra (see lWest et al.ll2004l . for more details). The fact that we 



successfully recovered the same active stars gives us confidence that FLI values are a useful 
measure in cases where the EW is ill-defined, and that the FDR method can be used to 
distinguish between the distributions. 



4. Flare Criteria 



We used two separate methods to define whether or not an object in our sample con- 
tained a flare. The first was to use FDR analysis to determine which spectra had emission 
in both Ha and H/3 that was strong enough to be considered a flare. The other method was 
to find objects whose Ha and H/3 lines were both in emission and showed large variations 
together with time. These two methods are not mutually exclusive; flares that have strong 
emission lines that vary over the course of several exposures will meet both criteria. 



4.1. Strong Emission Lines using FDR Analysis 

We used the DR5 composite spectra of the active stars (i.e. those that met the emission 
threshold criteria determined from §3.3[) as the null distribution for the flare FDR analysis. 

Although the objects that make up the null distribution (the composite spectra) are 
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the same objects that make up the test distribution (the individual exposures), the distri- 
butions are distinct in two significant ways. First, the signature of small flares that are only 
present in a subset of the exposures is diluted in the composite spectra. Second, outliers 
were removed from the null distribution by rejecting spectra that were not part of the con- 
tinuous distribution of FLI values. A few small or diluted flares may still be present in the 
null distribution. However, any flares in the null will simply increase the threshold value 
determined in the analysis and decrease the number of objects classified as flares. 

As in the active sample determination described in §3.31 we first normalized the null 
and test distributions, and then applied the FDR analysis separately to the HSN and MSN 
samples, for each M subtype. Again we found that the HSN sample analysis was limited to 
MO-3 types and the MSN sample to MO-7 types. The flare threshold values for both Ha and 
H/3 are given in Table 2, and are shown for the HSN sample in Figure El To be classified as 
a flare, we required that an individual spectrum had FLI values for both Ha and H/3 above 
the respective FDR flare threshold values. Thus, many of the spectra above the Flare line in 
Figure |6] were not found to be flares because they did not make the H/3 cut. This effectively 
removed cases of Ha activity that were not strong enough to be classified as flares. 



4.2. Variable Emission Line Strength 



We also identified flares by looking for variation in emission line strength among indi- 
vidual exposures of the same objec t. H a emission is kno wn to vary outside of flaring (see 



discussion in §|3]); iGizis et al.l (120021 ) and iLee et al.l (120101 ) found that "quiescent" Ha emis- 



sion often varies by as much as 30% on short timescales (minutes-hours). This may well be 
due to low-level flaring, but it is not possible to classify unequivocally as flare emission due 
to the relatively poor time resolution of their data (and ours). We therefore concentrated 
on variability that exceeds this level (representing larger, definite flares), and developed the 
following flare criteria based on both the Ha and H/3 FLI values: 



1. Mean FLI values > 3 over all exposures for an object for both Ha and H/3; in these 
strong emission line spectra, we identified flares if both the Ha and H/3 FLI values 
for an object differed by more than 30% between the minimum and maximum values 
measured over the course of the exposures. 

2. Mean FLI values < 3 in both Ha and H/3; we identified flares in these weak emission 
line data using the criterion that max FLI - min FLI > 3. We adopted an absolute 
change criterion because a 30% change is still within the noise when the mean FLI 
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value is small. Since FLI is denned as I /a, this is essentially a 3a requirement (max 
FLI is 3a larger than min FLI). 

3. Mean FLI values > 3 in Ha but < 3 in H/3; these objects have moderate emission 
strength in Ha but weak emission in H/3. We require that they meet both of the above 
variability criteria (max FLI - min FLI > 3; percent change > 30%) to be considered 
a flare. 

We examined objects that narrowly missed passing these flare criteria and confirmed 
that these quantitative measures are robust, and no real flares were omitted from the sample. 
Examples of flares that were identified based on each of these three conditions are shown in 
Figure UJ A large, decay-phase flare (top left) that met the first variability criterion (VI) 
also met the FDR strong lines criterion. The example flare that met the second variability 
criterion (V2, top right) demonstrated weak emission during the first two exposures, but 
had significantly stronger emission in the final exposure. The moderate quiescent emission 
example (bottom left) met the third variability criterion (V3) with a flare during the third 
exposure. The flare that met the FDR strong lines criterion (bottom right) showed very 
little change between exposures. 



5. Results 

Our automatic flare identification using the above methods resulted in 285 individual 
spectra on 72 M dwarfs being classified as flaring. We visually examined each spectrum and 
rejected 9 objects due to calibration errors, white dwarf companions, or cataclysmic variable 
signature. Our final flare sample therefore comprises 243 spectra on 63 stars. The properties 
of the 63 flaring stars are listed in Table 3. 
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Table 2. FLI threshold values 
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Table 3. Properties of the flares in the SDSS Spectroscopic Sample 
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Table 3 — Continued 
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A summary of the number of stars and exposures that met each flare criterion is pre- 
sented in Table 4, which also shows that the majority of flares come from the MSN sample. 
This is not unexpected, since there are many more stars in that sample. Our results suggest 
that flares in this sparsely sampled spectroscopic data set are most easily identified through 
emission line variability rather than the absolute strength of the emission lines. 

5.1. Time Evolution of Flares 

We first analyzed our flare sample to verify that the properties of the flares we identified 
are consistent with the characteristic time evolution seen in continuous monitoring obser- 
vations. Because our data might capture any or all phases of a flare, we expected to see 
examples of flares in the rise, peak, and decay phases as well as a few large impulsive phase 
flares that show blue continuum flux enhancement. 



5.1.1. Rise, Peak, and Decay Phase Flares 

We defined decay-phase flares as those showing both Ha and H/3 in at least two consec- 
utive exposures with diminishing FLI values. Rise-phase flares were defined in the opposite 
sense, with two or more consecutive exposures having increasing FLI values for both Ha and 
H/3. There were ten rise phase flares and 11 decay phase flares in the sample. Additionally, 
we saw five "peak" flares that show first increasing and then decreasing FLI values, which we 
interpret as flares that were seen through most of their evolution. We also identified flares 
that occurred in only a subset of the exposures, either in the final exposure of a sequence, 
or in cases where some of the exposures are separated by many hours or days from the 
other exposures. In 19 cases, at least one exposure was obtained that showed a quiescent 
spectrum with no sign of flaring. Flares with quiescent spectrum, as well as rise, decay, and 
peak phase flares are noted in the final column of Table 3. The time evolution of the flares 
we have identified is consistent with the known characteristics of flares. 



5.1.2. The Balmer decrements of Rise and Decay Flares 

Previous observations using time resolved spectroscopic monitoring have shown that 
the higher order Balmer lines show larger increases in flux during flares than the lower order 
lines. This results in the observed Balmer decrement (ratio of individual lines to a fiducial, 
often H/3 as we adopt here) becoming flatter, since H7 and the higher order lines strengthen 
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Table 3 — Continued 



RA Dec 


SpT 


Sample 


# exp. a 


ID Method b 


Notes c 
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b The method used to identify flares. SL is the strong line criterion (see 
§4.ip and V1,V2, and V3 are the three variability criteria (see §4.21) 

c "Rise" , "peak" , and "decay" phase flares are defined in §5.1.11 We also 
note flares in which at least one spectrum is in quiescence, as well as the 
two flares with enhanced continuum emission (discussed in §5.1.3p . 



Table 4. Flare identification by emission line strength and variability criteria 
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relative to H/3 (IGarcfa- Alvarez et al.ll2002r lAUred et al.ll2006l see also Figure 7 of Bochanski 
et al. 2007) This behavior can be explained by a model where flare heating increases the 
local electron density at chromospheric temperatures. The higher order Balmer lines are 
less optically thick than the lower order lines, so the increase in electron density causes 
a relatively larger increase in their emission. The Ha line behavior relative to the other 
lines has previously been more difficult to measure d uring flares, since observat ions typically 
concentrated on the blue part of the s pectrum (e.g. lHawley fc Pettersenlll99ll ). Along with 
their blue spectra, (jEason et al.lll992l ) obtained data with a different instrument at lower 
time resolution showing that the Ha line, being alrea dy very optically th ick, does not respond 
strongly during a large flare in accord with models ( lAUred et al.l 120061 ). 



Figure E] illustrates the line flux normalized to the bolometric luminosity for several 
lines (top), as well as Balmer decrements for four flares (bottom; two rise- and two decay- 
phase). The Balmer lines increase during the rise phase flares, and decrease during the 
decay phase flares, as expected. In the bottom panels, we see that for all of the Balmer lines 
except Ha, the decrements are flat throughout the flares. As expected, the Ha decrement 
decreases during the rise-phase, when H/3 is becoming relatively stronger, and increases 
during the decay-phase, as the atmosphere gradually returns to its quiescent temperature and 
density structure. Our observations are in accord with expectations, and actually represent 
truly simultaneous measurements with the same instrument of the entire Balmer decrement 
(including Ha) during stellar flares. 

The Ca II K line r esponds less s trongly, and rises and deca ys more slowly during flares, 
than the Balmer lines (IButlerl Il991t lHawley fc Pettersenlll99ll ). as is also evident from our 
observations. 



5.1.3. Continuum Flux Enhancement 

Large flares cause strong blue continuum flux enhancement, with less dramatic effects 
at redder wavelengths. Only very strong flares caught in the impulsive phase would lead to 
a significant flux enhancement for the wavelength range and relatively long exposure times 
(which dilute the short-lived continuum burst) of the SDSS spectra. We expected to find only 
a small number of such flares, since large flares are rare and we are only sensitive to measuring 
continuum enhancement on HSN stars or very bright MSN stars, where we have a good 
measure of the quiescent continuum. Our sample contained two continuum flares, an example 
of which is shown in Figure |9j Both flares occurred during the third and final exposure of a 
series, and both show very strong emission lines. These are evidently large flares that were 
caught during the impulsive phase. We measured the flare-only spectrum by subtracting 
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the quiescent (first) exposure obtained for each star from the flare exposure. We found that 
the continuum flux represented 88% and 90% of the total flare radiation, respectivel y , with 
the remainder of the energy coming from the line radiation. lHawley fc Pettersenl (Il99ll ) 
found that the continuum emitted 96% of the flare flux during the impulsive phase of their 
very large flare, and averaged 91% over the course of the flare. This suggests that our flare 
exposures may have included the impulsive phase and part of the decay phase of each flare. 
The existence of these continuum flare spectra in our sample, and the relative strength of 
the lines and continua, give us additional confidence that we are correctly identifying flares 
by our automatic algorithms. 



5.2. Flare Duty Cycle 

Figure [10] shows our observed flare duty cycle as a function of spectral type. The duty 
cycle increases with later spectral type, from 0.02% for M0 stars to 20% for M9 stars. The 
error bars were calculated from binomial statistics, and the number of flares identified at 
each spectral type is noted. The numbers in parentheses are the number of flares identified 
by the variability criteria. The much larger duty cycles for M7 and later are based on a 
small number of stars, and although there are many thousands of objects in the M0-M3 
bins, there are very few flares. In order to improve the statistics, we binned by spectral type 
and found the duty cycles to be 0.022 ±0.016% of M0-1 stars; 0.077 ±0.023% for M2-3 stars; 
0.56 ± 0.09% for M4-6 stars; and 3 ± 1% of M7-9 stars. 



Kowalski et al.l (120091 ) investigated flares in a large SDSS photometric sample. They 
required near simultaneous photometric enhancements in both the u and g filters, and em- 
ployed a threshold of \Au\ > 0.7 to define a flare. Figure ITT1 compares their flare duty cycle 
with our results for the binned data. We find a higher duty cycle, by a factor of ~ 4.5, 
but the duty cycle as a function of spectral type has a similar trend in both studies. It 
is not surprising that the duty cycles are different, since the two studies use such different 
observations and threshold criteria. However, we can qualitatively estimate the effect of 
observing line versus continuum enhancements, and the effect of exposure time. The line 
emission remains enhanced for much longer than the continuum emission leading to a higher 
probability of observing flares in line emission. Longer exposure times increase the chance 
of a flare occurring during the exposure. Both of these factors lead to the observed higher 
duty cycle seen in our spectroscopic observations. 
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5.3. The Vertical Distribution of Flaring Stars in the Galaxy 



Nearly all fl a res occur on magnet ically active stars, as defined by Ha in emission 
(IReid et al.l Il999l ; iKowalski et al.l 120091 ) . As discussed in Section 13.31 the active fraction 



increases with spectral subtype. In addition, IWest et al.l (120081 ) found that the active frac- 
tion decreases with vertical height from the Galactic Plane (denoted by |Z|) for all spectral 
types, which they attribute to an age effect, where older stars are less likely to be active. 
The age at which activity disappears changes with spectral type however, with later type M 
dwarfs remaining active for a longer time. 

We limited our analysis of the vertical distribution of the flaring stars to just the M4-6 
dwarfs, which show a large number of flares, and have a similar flaring fraction (see Figure 
ITUl) . Earlier type stars have very few flares, and later type stars are too faint to be observed 
at large distances. 



Figure [T21 (top panel) shows the distribution of active M4-6 dwarfs and the fraction that 
showed flares as a function of vertical distance from the Galactic Plane. We repeated our 
analysis, adjusting both the bin size and bin positions and verified that the increase in flare 
fraction as a function of |Z| is a real effect and not an artifact of the binning process. The 
bottom panel of Figure [T21 shows the cumulative distributions of the flare stars (black line), 
the active stars (dotted gray line) and the inactive stars (dashed lin e). Although th e re are 
thousands of obje cts be yond | Z | =200 pc, none of them showed flares. IKowalski et al.l (120091 ) 
(120071 ) also found that all of the flares on M4 dwarfs occur within 200pc of 



and Welsh et al 



the Plane. 

Our results indicate that flares occur on stars that are preferentially closer to the Plane 
than the inactive stars, and closer even than the active population. If we interpret this in 
terms of an age effect, it suggests that flare stars (which are a subset of the active population) 
are younger than the magnetically active stars as a whole. Evidently the frequency of flaring 
decreases as a star ages, although it still may exhibit quiescent chromospheric emission and 
therefore be counted as active. 



6. Conclusion 

We have presented the first study to systematically search for flares in a sample of 
tens of thousands of time-resolved spectroscopic observations from the SDSS Data Release 
6. We developed two new techniques for identifying and analyzing flares in relatively low 
SNR spectra: the Flare Line Index (FLI) measurement and the application of the False 
Discovery Rate to flare spectra. FLI values provide a measure of emission line strength even 
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when the continuum level is not well-measured. The FDR analysis distinguishes between 
the overlapping distributions in line strength of strong quiescent emission and genuine flares, 
providing a threshold FLI value for being considered a flare, while controlling the rate of 
false-positives. Using our measured FLI values for the Ha and H/3 lines, we applied the FDR 
analysis and a set of variability criteria to find flares. 

We identified flares on 63 stars from the 30,971 total stars comprising both the HSN 
and MSN samples. The flares include examples of rise and decay phases, as well as some 
sequences of exposures that appear to trace the flare through its entire evolution. The 
Balmer decrements are consistent with results from time resolved spectroscopic monitoring 
of individual flares and expectations of models. We also found two flares that show continuum 
flux enhancements, which occur during the initial impulsive phase of large flares. Because the 
properties of the flares in our sample are consistent with previous findings, we are confident 
that our techniques, which we developed specifically to work with coarsely time-resolved 
spectroscopic observations, successfully identified flares. 

We found that the flare duty cycle increases monotonically as a function of M dwarf 
spectral subtype, from 0.022 ± 0.016% of M0-1 stars to 3 ± 1% of M7-9 stars. The flare 
distribution as a function of spectral subtype is similar to the distribution of active stars, 
providing further evidence that most flares occur on active stars. Nearly all of the flares in the 
sam ple occur on stars nea r the Galactic Plane, a result that is consistent with other studies 



stars closer to the Plane than the mean active population, indicating that flare frequency 
likely decreases with age. 

This study represents an important step toward understanding the intrinsic flare fre- 
quency distribution as a function of stellar active lifetimes. Surveys that don't have high 
time resolution and continuous coverage (from which the total flare energy can be computed) 
require detailed and nuanced Monte Carlo simulations in order to constrain the physical flare 
frequency distribution. Simulations which consider stellar samples in the Galactic context is 
the subject of ongoing study (Hilton et al., in prep). 
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Fig. 1. — Histograms of the exposure times for all component spectra in the sample (left) 
and of the exposure times fo r the component spectra (right). The dashed histograms are for 
the entire IWest et al.l ( 120081 ) sample, while the solid histograms are the stars considered in 
this study. 
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Fig. 2. — A sequence of spectra of a typical decay phase flare (left), with the time in 
minutes from the beginning of the first exposure. A detailed look at the emission lines 
(right, color corresponding to the exposure number on the left) reveals that Ha and H/3 
remain in emission longer than the higher-order Balmer lines. Note that Ca II K increases 
during the third exposure, an effect not observed in the Balmer lines. 
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Fig. 3. — The distribution of spectral subtypes for the HSN (dotted), MSN (dashed), and 
total (solid) samples. As expected, the HSN sample mostly contains the earlier, and intrinsi- 
cally brighter, spectral subtypes. Although the MSN sample has more later spectral subtype 
stars, there are relatively few beyond M6. 
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Fig. 4. — Three consecutive spectra (ordered by time: black/solid, red/dotted, blue/dashed) 
of an M5 star from the MSN sample that show no evidence of strong line flux in the H/3 
(top) and Ha (bottom) line regions. The vertical lines mark the line regions used in the 
calculations. The EW values for the H/3 are spuriously high in some exposures, while all of the 
FLI values are low. For the Ha line, the FLI values are similar to the EW measurements. This 
demonstrates the usefulness of using FLI values when the continuum is weak or uncertain. 
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Fig. 5. — The distribution of FLI values for Ha of the component spectra of M3 stars in the 
HSN sample (gray line) and the distribution of the null sample taken from the inactive DR5 
integrated spectra (black line). The FDR-derived threshold value (vertical line) separates 
the inactive from the active and flaring spectra, which are seen as the high FLI value tail of 
the HSN sample. 
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Fig. 6. — The Ha FLI values from the HSN sample as a function of spectral subtype are 
shown as gray points. Each point has been given a random offset within its spectral type bin 
and enlarged where the point density is low for clarity. Most spectra have low FLI values, 
with the mean increa s ing w i th la t er spe ctral type, as expected from the activity fraction 
studies of IWest et all fcooi bood . l2008h . The solid line labelled "Emission" indicates the 
FDR-derived threshold values separating active spectra from inactive ones, while the "Flare" 
line denotes the threshold values separating active stars from flares. Because our flare criteria 
require both Ha and H/3 to be above the relevant threshold value, not all points above the 
"Flare" line are categorized as flares in our sample. 
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Fig. 7. — The H/3 and Ha lines of typical flares that meet each of the flare criteria (see §4.11 
and 14.21 for details about the strong lines and variability criteria). Each panel shows all of 
the exposures for that object, with the black, red, and blue lines corresponding to the first, 
second, and third exposures, respectively. Note that the flare that meets the first variability 
criterion also meets the FDR strong lines criterion. 
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Fig. 8. — Line fluxes of emission lines normalized by the bolometric luminosity for two rise 
phase flares and two decay phase flares (top). Note that Ca II K (crosses) does not follow 
the same trend as the Balmer lines. The bottom panels show the Balmer decrements for the 
same flares. The Ha decrement is smaller when the flare is closer to its peak (later for rise 
phase flares, and earlier for decay phase flares). 
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Fig. 9. — The spectrum of an M6 dwarf (R.A., Dec. = 194.47690, 3.55015) during quiescence 
(black) and during a flare (blue). This exposure captures at least part of the impulsive phase 
of the flare and shows a clear continuum enhancement in the blue portion of the spectrum 
(shown in the bottom panel), as well as strong emission lines. 
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Fig. 10. — The flare duty cycle as a function of spectral type. The error bars are calculated 
from the binomial distribution. The number of flares in each bin is given, with the number 
of flares that meet the variability criteria given in parentheses. The flare fraction increases 
with spectral subtype as expected from previous studies, and is similar to the rise in the 
active fraction of M dwarfs. 
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Fig. 11. — The flare duty cycle as measured by lKowalski et al.l (120091 closed circles) and this 
study (stars), for spectral type bins M0-1, M2-3, M4-6, M7-9. Scaling the iKowalski et al. 
(120091 ) results by a factor of 4.5 (squares) makes them consistent with our results. 
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Fig. 12. — Top: The flare duty cycle as a function of dista nce from the Gal actic midplane 
(|Z|) for M4-6 stars (black line). The number of active stars (IWest et al.ll2008l ) in the sample 
is also shown (gray line, dotted, right axis). The fraction of objects that flare decreases 
rapidly with distance from the Galactic Plane. Bottom: the cumulative distribution of 
absolute Z distance for the flare, active, and inactive stars between subtype M4 and M6. 
These results indicate that flaring occurs preferentially closer to the Galactic Plane than the 
general active star population, and therefore on younger stars. 



